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ABSTRACT 


A  "hydrogen  partitioning"  model  has  been  developed  to  account  for  the 
pressure  and  temperature  dependence  for  hydrogen-assisted  crack  growth.  The 
model  gives  explicit  recognition  to  the  role  of  hydrogen-microstructure  inter¬ 
actions  in  determining  the  distribution  (or  partitioning)  of  hydrogen  among 
the  various  microstructural  elements  (principally  between  the  prior-austenite 
grain  boundaries  and  the  matrix)  and  the  rate  of  crack  growth  along  the 
elements.  It  also  takes  Into  account  the  role  of  various  rate^wntrol  1  Ing 
processes  In  determining  the  rate  that  hydrogen  is  being  suppl  lett-tQthe 
fracture  process  (or  embrittlement)  zone. 

Quantitative  assessment  of  the  model  Indicates  very  good  agreements 
between  the  model  predictions  and  the  observed  crack  growth  responses  for  AISI 
4340  and  4130  steels  tested  in  hydrogen  and  for  AISI  4340  steel  tested  in 
hydrogen  sulfide.  This  model  accurately  characterizes  the  reduction  in  crack 
growth  rate  and  the  concomitant  change  In  fracture  mode  at  "high*  tempera¬ 
tures.  Through  its  Integration  with  the  earlier  models,  based  on  rate  con¬ 
trolling  processes,  the  model  predicts  the  pressure  and  temperature  dependence 
for  K- Independent  crack  growth  over  the  entire  range  of  environmental  condi¬ 
tions. 
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1.  INTRODUCTION 

Subcritical  crack  growth  in  high  strength  steels  exposed  to  gaseous 
environments  has  been  widely  studied  from  an  interdisciplinary  point  of 
view.  A  series  of  parallel  fracture  mechanics  and  surface  chemistry 
studies  on  AISI  4130  and  4340  steels,  exposed  to  hydrogen  containing 
gases  (such  as,  hydrogen,  hydrogen  sulfide  and  water  vapor)  and  to  water, 
have  provided  a  fair  understanding  of  environmentally  assisted  crack 
growth  [1-11]. 

It  is  now  clear  that  environmentally  assisted  crack  growth  is  con¬ 
trolled  by  a  number  of  processes  in  the  embrittlement  sequence  [4].  The 
various  processes  that  might  be  involved  are  (i)  transport  of  the  gas  or 
gases  to  the  crack  tip,  (ii)  the  reactions  of  the  gas  or  gases' with  newly 
created  crack  surfaces  to  evolve  hydrogen  (viz.,  physical  and  dissociative 
chemical  absorption  in  sequence),  (ill)  hydrogen  entry  (or  absorption), 

(iv)  diffusion  of  hydrogen  to  the  fracture  (or  embrittlement)  sites,  and 

(v)  hydrogen-metal  interactions  leading  to  embrittlement  (i.e.,  the 
embrittlement  reaction),  see  Fig.  1.  According  to  the  multi -reaction 
rate  theory  [12],  the  overall  crack  growth  response  is  governed  by  the 
slowest  unit  process  in  the  chain  of  processes  operating  in  sequence. 

Any  one  of  these  processes  may  act  as  the  rate-controlling  process, 
depending  on  the  environment-material  condition  (material,  environment, 
pressure,  temperature,  etc.). 

For  AISI  4130  and  4340  steels  exposed  to  hydrogen  [1-3,  6-8]  and 
hydrogen  sulfide  [9],  the  K-independent  Stage  II  crack  growth  exhibits 
substantially  different  responses  In  two  temperature  regions.  Fig.  2. 

At  "low"  temperatures ,  crack  growth  reflects  control  by  gas  transport  to 


the  crack  tip  (hydrogen  sulfide  at  the  lower  pressure,  curve  b),  by  sur¬ 
face  reaction  (hydrogen  or  water,  curves  c  and  d),  and  by  hydrogen  dif¬ 
fusion  (hydrogen  sulfide  at  the  higher  pressure,  curve  a).  The  crack 
growth  rates  conform  to  the  temperature  and  pressure  dependence  of  the 
rate  controlling  process  [11].  At  "high  temperatures",  the  growth  rates 
are  substantially  lower  than  those  given  by  these  processes,  and  some  other 
process  appears  to  be  in  control.  The  intermediate  region  (Region  B) 
represents  the  transition  between  Regions  A  and  C.  The  fundamental  ques¬ 
tion  posed  by  these  results  Is  which  process  or  processes  cause  the  de¬ 
crease  in  crack  growth  rate  at  "high  temperatures". 

Attempts  to  explain  this  behavior  were  first  made  by  Williams, 

Nelson  and  Tetelman  [1-3]  in  the  early  1970's.  According  to  their  surface 
reaction  model,  the  high  temperature  region  of  crack  growth  response  is 
controlled  by  the  equilibrium  between  absorbed  hydrogen  and  gaseous  hydro¬ 
gen.  Subsequently,  several  "surface  reaction"  models  were  proposed  by 
other  researchers  [13-15].  These  models  are  purely  chemical  and  considered 
only  surface  reactions  in  controlling  crack  growth  response  even  though 
other  processes  might  have  been  In  control  (Fig.  2)  [9-11], 

Recent  fractographic  studies  on  an  AISI  4340  steel  indicated  that 
the  reduction  In  crack  growth  rate  with  Increasing  temperature  into  the 
high  temperature  region  is  accompanied  by  a  change  in  fracture  surface 
morphology;  namely,  from  Intergranular  to  transgranular  (micro-void 
coalescence  type)  failure  [16].  This  finding  implied  that  the  associated 
change  In  hydrogen  supply  (as  predicted  by  the  surface  reaction  models) 
could  not  be  the  sole  cause  for  the  observed  change  in  crack  growth 
response.  The  observed  difference  In  fracture  surface  morphology 


suggested  that  different  micromechanisms  or  different  hydrogen  embrittle¬ 
ment  processes  might  be  in  control  for  the  two  temperature  regions. 

To  clarify  this  issue,  the  response  of  different  microstructural 
elements  to  hydrogen  embrittlement  and  the  partitioning  of  hydrogen  among 
these  elements  as  a  function  of  temperature  and  pressure  are  considered. 

A  "hydrogen  partitioning"  model  for  crack  growth  in  gaseous  environments 
is  proposed.  The  model  provides  a  direct  correlation  of  crack  growth 
kinetics  with  the  response  of  microstructural  elements  to  hydrogen  embrittle¬ 
ment. 

2.  EXPERIMENTAL  BASIS  FOR  THE  MODEL 

The  proposed  model  is  based  on  the  following  experimental  findings 

[16]: 

(1)  In  the  low  temperature  region,  the  fracture  surface 

morphology  is  characterized  by  intergranular  separation 
along  prior-austenite  grain  boundaries,  with  a  small 
amount  of  transgranular  quasi -cleavage  [2,3,17,18]. 

The  fracture  mode  is  unaffected  by  the  rate  controlling 
process  for  crack  growth  as  long  as  an  adequate  amount 
of  hydrogen  is  supplied. 

(11)  The  change  in  Stage  II  crack  growth  rate  in  going 
from  low  to  high  temperature  ,  with  a  maximum 
rate  at  an  Intermediate  temperature,  is 
related  to  a  transition  from  intergranular  and  quasi- 
cleavage  to  microvoid  coalescence  modes  of  separation. 


(iii)  The  temperature  at  which  the  transition  in  fracture 
surface  morphology  occurs  is  dependent  on  hydrogen 
and  hydrogen  sulfide  pressure.  The  higher  the 
gas  pressure,  the  higher  is  this  temperature.  This 
trend  is  consistent  with  the  effect  of  gas  pressure 
on  the  change  in  crack  growth  response  from  low  to 
high  temperatures. 

These  fractographic  observations  are  difficult  to  explain  in  terms  of 
the  models  proposed  by  Williams  and  Nelson  et  al.  [1-3,  13-15].  The 
observed  difference  in  fracture  surface  morphology  and  the  associated  reduc¬ 
tion  in  crack  growth  rates  suggest  that  different  micromechanisms  or  dif¬ 
ferent  hydrogen  embrittlement  processes  need  to  be  considered  for  the  two 
temperature  regions. 

To  examine  the  role  of  micromechanism,  one  might  divide  the  embrittle¬ 
ment  sequence  roughly  into  three  portions  [16].  These  three  portions  are 
depicted  schematically  in  Fig.  3.  The  first  portion  is  concerned  with  the 
supply  of  hydrogen,  and  includes  gas-phase  transport  (i.e.,  external  trans¬ 
port),  surface  reactions,  and  diffusion  (i.e.,  internal  transport)  of  hydro¬ 
gen  to  the  embrittlement  sites.  The  second  portion  involves  the  partition¬ 
ing  of  hydrogen  among  the  potential  fracture  sites  in  the  microstructure. 
These  sites,  in  AISI  4340  steel,  include:  (i)  the  prior-austenite  grain 
boundaries,  (ii)  the  martensite  lath  or  patch  boundaries,  and  {110}  ,  and 

ct 

{112>ai  planes  through  the  martensites,  and  (iii)  the  martensite  lattice 
itself.  The  third  portion  deals  with  the  embrittlement  reaction  at  each 
type  of  sites,  with  the  prior-austenite  grain  boundaries  being  most  suscep¬ 
tible  and  the  martensite  lattice  least  susceptible  to  hydrogen  embrittlement 


Cracking  at  the  respective  mi crostructural  elements  would  depend  on  the 
concentration  and  rate  of  supply  of  hydrogen  there.  The  overall  crack 
growth  rate  would  be  determined  by  the  rates  of  cracking  through  the  par¬ 
ticipating  mi crostructural  elements. 

It  is  reasonable  to  suggest  that  the  partitioning  or  distribution  of 
hydrogen  to  the  different  mi crostructural  elements  is  a  function  of  temper¬ 
ature  and  pressure.  Hence,  the  concentration  and  the  rate  of  hydrogen 
supplied  to  the  different  fracture  sites  become  functions  of  temperature 
and  pressure.  At  low  temperature  (Region  A),  the  prior-austenite  grain 
boundaries  and  the  slip  planes  are  saturated  by  hydrogen.  The  fracture  is 
therefore  predominantly  intergranular,  and  includes  a  small  amount  of 
quasi -cleavage.  Because  the  rate  of  supply  of  hydrogen  from  the  external 
environment  is  relatively  slow,  crack  growth  is  controlled  by  one  of  the 
hydrogen  supply  processes  [10,11].  At  the  higher  temperatures,  or  Region 
C,  the  hydrogen  concentration  at  grain  boundaries  and  slip  planes  is  expec¬ 
ted  to  drop  quickly  and  more  hydrogen  goes  into  the  martensite  lattice. 
Increasing  amounts  of  microvoid  coalescence  or  dimpled  failure  now  occur 
with  increased  temperature  and  result  in  slower  crack  growth  rates. 
Accordingly,  the  changes  in  crack  growth  rate  and  crack  growth  response 
with  temperature  are  attributed  now  to  the  transfer  of  micromechanisms  of 
fracture,  instead  of  the  processes  of  hydrogen  supply.  These  considerations 
are  equally  applicable  to  hydrogen  and  hydrogen  sulfide. 


3.  ASSUMPTIONS  AND  BASIC  EQUATIONS 

In  developing  the  model,  the  following  five  assumptions  are 

made: 

{ i )  Crack  growth  is  at  steady  state,  and  is  in  Stage  II.  Hence, 
the  crack  growth  rate  remains  constant  at  a  given  temper¬ 
ature  and  hydrogen  pressure,  and  is  independent  of  the 
stress  intensity  factor  K. 

(ii)  The  rate  of  hydrogen  supplied  to  the  fracture  process 
zone  is  governed  by  the  controlling  process  for  hydro¬ 
gen  supply.  These  processes  Include  (a)  transport  of 
the  deleterious  gas  to  the  crack  tip,  (b)  surface  reaction 
of  the  gas  with  the  newly  created  surface  at  the  crack 
tip,  and  (c)  diffusion  of  hydrogen  to  the  embrittlement 
region  ahead  of  the  crack  tip. 

(ill)  The  controlling  process  of  hydrogen  supply  remains 
unchanged  over  the  range  of  temperature  of  interest; 
that  Is,  from  Region  A  to  Region  C. 

(1v)  The  rate  of  crack  growth  for  each  fracture  mode  is 
determined  by  the  rate  of  supply  of  hydrogen  to  each 
type  of  fracture  site,  Q^. 

(v)  The  different  micromechanisms  of  cracking  operate  in 
parallel  or  concurrently. 

Based  on  these  assunptions,  the  crack  growth  rate  for  the  ith  fracture 
mode  (or  micromechanism),  (da/dt)^,  is  given  by  Eqn.  (1).  In 
the  equation,  a*  is  the  proportionality  constant  for  the  ith  mode. 


For  concurrent  growth,  the  mean  crack  growth  rate,  (da/dt)jj  is  the  sum 
of  the  products  of  growth  rate  and  the  corresponding  areal  fraction,  f j , 
of  the  individual  modes  on  the  fracture  surface,  and  is  given  by  Eqn.  (2). 

(da/dt)n  =  Ef1aiQi  (2) 

For  AISI  4340  steel  fractured  in  hydrogen  and  hydrogen  sulfide,  the 
areal  fraction  of  quasi-cleavage  (f^)  is  small  (less  than  about  14%)  as 
compared  with  that  of  grain  boundary  separation  at  low  temperatures, 
and  is  nearly  zero  at  the  higher  temperatures  [16].  The  contribution  of 
quasi -cleavage  to  the  main  crack  growth  rate,  therefore,  may  be  neglected. 
Then,  Eqn.  (2)  can  be  simplified  as  follows: 

(da/dt)ji  =  fb<»bQb  +  VjQ, 


*  <fbVb  *  fiVtw 


(3) 


where  Q  is  the  total  rate  of  hydrogen  supplied  to  fracture  process  zone, 
and  *cb,  are  the  hydrogen  distribution  coefficients  (<b  *  Qb/Q»  = 
Qj/Q,  and  *  1). 

Equations  (2)  and  (3)  are  two  of  the  basic  equations  for  the  model. 

For  simplicity,  only  Eqn.  (3)  is  used  in  further  quantification  of  the  model 


4.  DISTRIBUTION  COEFFICIENTS  AND  HYDROGEN  SUPPLY  RATE 


4.1  Distribution  Coefficients  and  k 

The  distribution  coefficients  tcb  and  are  defined  as  the  fraction 
of  the  rate  of  hydrogen  being  supplied  to  the  grain  boundaries  and  to  the 
lattice  respectively,  i . e. : 


*  VQ  •  (VV/'1  +  vv 

(4. a) 

•  qt/Q  -  1/(1  +  Qb«*> 

(4.b) 

To  estimate  these  coefficients,  hydrogen  distribution  in  a  moving 
frame  of  reference,  with  the  origin  at  the  crack  tip,  is  considered  first. 
Because  crack  growth  is  assumed  to  be  at  steady  state,  with  constant  speed 
(da/dt) j j ,  the. average  concentration  of  hydrogen  In  the  fracture  process 
zone,  Cy,  is  independent  of  time  within  this  moving  frame.  The  average 
hydrogen  concentration  in  the  grain  boundaries  (Cb)  and  in  the  lattice 
(C^)  in  the  same  zone  are  also  time  independent.  The  exact  location, 
shape  and  size  of  the  fracture  process  zone,  however,  cannot  be  specified 
at  the  present  time.  For  these  estimates.  It  is  sufficient  to  assume  it 
to  be  inside  a  circle  of  diameter  with  its  center  at  the  crack  tip. 

Fig.  4,  and  to  assume  that  hydrogen  is  confined  within  this  circle. 

It  can  be  readily  shown  that  the  rates  of  hydrogen  supply  to  the 
fracture  process  zone,  and  to  the  grain  boundaries  and  the  lattice  within 
this  zone  are  given  by  Eqns.  (5a)  to  (5c),  respectively: 

Q  =  CvB*(da/dt)n  (5. a) 


Qh  -  Ch$Bfc(da/dt)TT 


(5.b) 


Q4  =  Ca(l-6)B*(da/dt)n 


(5.c) 


In  Eqn.  (5),  6  and  (1-6)  are  the  volume  fractions  of  prior-austenite 
grain  boundaries  and  lattice,  respectively  in  the  fracture  process  zone. 
Hydrogen  concentration  is  given  in  moles  per  unit  volume.  By  substituting 
(5.b)  and  (5.c)  into  (4. a)  and  (4.b),  the  distribution  coefficients  become 


C«/(i-«)](cb/c£) 


(6. a) 


(6.b) 


In  general,  6  is  negligibly  small  compared  with  1.  Therefore,  Eqn.  (6. a) 
and  (6. b)  can  be  further  simplified  to  the  following  form: 


Kb  i+6y^ 


(7.a) 


'  UICb/Ci 


(7.b) 


The  distribution  coefficients  and  k^)  are  now  given  in  terms  of  the 
ratio  of  local  concentrations  of  hydrogen  in  the  grain  boundaries  and 
in  the  lattice,  and  the  volume  fraction  of  grain  boundaries. 

It  is  well  demonstrated  that  grain  boundaries  are  strong  trapping 
sites  for  hydrogen  [19-21].  The  equilibrium  partitioning  of  hydrogen 
between  the  grain  boundaries  and  the  lattice  can  be  expressed  by  either 
Boltzmann  statistics  (for  dilute  solution  [22])  or  by  Fermi-Dirac  statistics 
(for  more  concentrated  solution,  or  at  low  temperatures  [23]).  Even  though 
crack  growth  is  assumed  to  be  at  steady  state,  conditions  of  equilibrium 


may  not  always  be  established.  Following  Bockris  [24],  therefore,  a  "non- 
equilibrium"  parameter  t  is  incorporated  in  the  Boltzmann  and  the  Fermi  - 
Dirac  statistics: 

=  t(a3/n)  Nx  exp  (Hg/RT)  Boltzmann  (8. a) 

t(a3/n)Nx  exp (Hg/RT) 

Cb/C*  °  r+ii(aVn)«Afil  exp(yftT)  Fermi-Dirac  (8-b) 

The  parameter  t  is  equal  to  1  when  trapping  is  in  equilibrium.  Nx  is  the 
density  of  trap  sites  in  the  grain  boundaries  (number  of  sites/m  ),  a  is 
the  lattice  parameter  of  the  metal,  n  is  the  number  of  atoms  per  unit  cell, 
is  Avogadro  constant,  and  Hg  is  the  binding  enthalpy  of  hydrogen  to  the 
grain  boundary.  Detailed  derivation  of  Eqns.  (8. a)  and  (8.b)  is  given  in 
Appendix  I.  Final  expressions  for  and  can  be  readily  obtained  by 
substituting  Eqns.  (8. a)  and  (8.b)  into  Eqns.  (7. a)  and  (7 . b) . 

It  should  be  noted  that  either  the  Boltzmann  or  the  Fermi -Dirac 
statistics  may  be  used.  The  choice  between  these  two  statistics,  however, 
is  not  crucial  to  modeling  and  quantitative  analyses,  because  both  sets 
give  the  same  values  for  and  <z  at  low  temperatures*.  At  high  tempera¬ 
tures,  Boltzmann  and  Fermi-Dirac  statistics  themselves  are  essentially  the 
same.  For  simplicity,  Boltzmann  statistics  are  chosen  for  use  in  the 
present  study. 

*  — 

By  examining  Eqns.  (8. a)  and  (8.b),  it  can  be  seen  that  Cb/Ce»l  when 
the  temperature  is  low  for  both  statistics.  Substituting  k  into  Eqns. 

(7. a)  and  (7.b)  gives  <b3l,  and  a  negligibly  small  value  for 


•  *•  > 
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4.2  Hydrogen  Supply  Rate  Q 
According  to  assumption  (2),  the  rate  of  hydrogen  supply  is  controlled 
by  one  of  the  processes  for  hydrogen  supply;  with  a  particular  process  in 
control  over  a  given  range  of  pressure  and  temperature.  Using  results 
from  recent  studies  and  from  modeling  of  crack  growth  for  each  of  the  con¬ 
trolling  processes  [4-11,  25],  the  hydrogen  supply  rate  (Q),  for  simplified 
cases,  may  be  given  in  the  following  form: 

For  transport  control  (based  on  Knudsen  flow)  [9,26] 

0  -  P0/T*  (9. a) 

For  surface  reaction  control  [8,  27] 

Q  -  exp(-Es/RT)  (9.b) 

For  diffusion  control  [28] 

Q  «  p0*  exp(-Ed/2RT)  (9.c) 

In  Eqn.  (9),  pQ  is  the  external  pressure;  m  is  the  exponent  for  pressure 
dependence  for  surface  reaction;  E$  and  Ed  are  the  activation  energies 
for  surface  reaction  and  for  diffusion,  respectively;  T  Is  the  absolute 
temperature;  and  R  is  the  gas  constant.  The  precise  form  of  these  equa¬ 
tions  depends  on  the  details  of  the  various  processes. 

5.  FORMULATION  OF  STAGE  II  CRACK  GROWTH  RATE 

By  incorporating  the  various  parameters  into  Eqn.  (3),  the  final 
formulas  for  Stage  II  crack  growth  are  obtained. 


6.  COMPARISON  OF  THE  MODEL  TO  EXPERIMENTAL  DATA 

Data  on  AISI  4340  steel,  tested  In  hydrogen  and  in  hydrogen  sulfide 
[7,9],  are  used  to  test  the  model.  Results  from  the  two  environments 
are  considered  separately  in  the  following  subsections. 

6.1  Hydrogen 

For  tests  in  hydrogen,  the  rate  of  hydrogen  supply  is  controlled  by 
the  rate  of  surface  reaction  [7].  In  this  case,  the  crack  growth  rate 
Is  given  by  Eqn.  (10. b).  Nearly  all  of  the  parameters  in  this  equation 
are  either  known  or  can  be  estimated  from  the  experimental  data.  Evalua¬ 
tions  of  these  parameters  are  described  in  Appendix  II. 

The  two  remaining  unknowns  are  the  density  of  trap  sites  In  the 
grain  boundaries  (Nx)  and  the  non-equilibrium  parameter  (t).  For  a  given 

19  23  _3 

microstructure,  Nx  must  be  a  constant  and  ranges  from  10  to  10  m 
[23].  The  parameter  x,  which  was  introduced  into  the  Boltzmann  equation 
to  account  for  non-equilibrium,  may  or  may  not  be  a  constant.  If  equili¬ 
brium  is  attained,  then  x  would  be  constant  and  Is  equal  to  one.  In 
general,  however,  x  is  expected  to  be  a  function  of  temperature  and 
pressure.  For  now,  as  an  initial  effort,  x  is  assumed  to  remain  constant 
for  a  given  pressure.  The  pressure  and  temperature  dependence  for  x 
will  be  considered  further  In  Section  7. 

For  the  convenience  of  numerical  analyses,  an  "apparent"  density 

of  occupied  trap  sites  in  the  grain  boundaries  (Napp)  is  defined,  with 
Na__  *  tN  .  Since  all  of  the  parameters,  except  Nann  in  Eqn.  (10. b)  are 

QpP  a  app 

known,  a  non-linear  least-squares  curve  fitting  program  can  be  used  to 
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search  for  the  value  of  Napp  that  provides  the  best  fit  to  the  experimen¬ 
tal  data.  The  best  values  of  Napp  obtained  are  3.8  x  1022,  2.88  x  1022 
and  8.56  x  1021  m" 3  for  hydrogen  pressures  of  133,  57  and  13.3  kPa, 
respectively.  Based  on  these  values,  corresponding  theoretical  curves 
of  crack  growth  over  the  temperature  range,  extending  from  Region  A  to  C, 
are  computed  from  Eqn.  (lO.b)  and  are  shown  in  Fig.  5.  These  results 
indicate  that  the  theoretical  curves  are  in  good  agreement  with  the  experi¬ 
ment.  The  parameters  used  in  the  equation  are  listed  in  Table  1.  It  is 
interesting  to  note  that  the  best  values  of  Napp  are  appreciably  differ¬ 
ent  for  the  three  different  pressures.  The  "apparent"  density  of  occupied 
trap  sites  in  the  grain  boundaries  increased  with  increasing  hydrogen 
pressure.  Because  Nx  is  expected  to  remain  constant,  this  fact  clearly 
Indicates  that  t  Is  not  equal  to  1  and  that  non-equilibrium  condition  needs 

to  be  considered.  The  physical  meaning  of  t  and  N.nn  is  discussed  in 

app 

Section  7. 

6.2  Hydrogen  Sulfide 

In  the  case  of  hydrogen  sulfide,  the  hydrogen  supply  rate  is  con¬ 
trolled  by  gas  transport  at  low  pressure  (0.133  kPa)  [9].  At  a  higher  pres¬ 
sure  (2.66  kPa),  the  supply  rate  Is  controlled  by  gas  transport  at 
intermediate  temperatures  and  by  diffusion  at  the  lower  temperatures 
[11,  16].  Hence,  in  examining  the  transition  from  Region  A  to  Region  C 
growth,  it  is  sufficient  to  consider  only  the  gas  transport  controlled 
case  for  both  the  low  and  high  pressures.  Thus,  Eqn.  (10. a)  is  used  to 
predict  the  crack  growth  response. 


The  same  values  of  Hg,  a,  n,  and  5,  used  for  the  hydrogen  case,  are 
used  here,  Table  2.  The  same  procedures  were  used  to  determine  the  coeffi¬ 
cients  and  the  other  parameters.  The  apparent  densities  of  occupied 

trap  sites,  N-__  =  xNw,  are  estimated  to  be  3.6  x  1022  m"3  and  4.2  x  1023  m 
app  x 

for  H2S  pressures  of  0.133  kPa  and  2.66  kPa,  respectively.  Again,  the 
predicted  curves  of  crack  growth  response  are  in  good  agreement  with  the 
experimental  data,  as  shown  in  Fig.  6. 

7.  FURTHER  CONSIDERATION  OF  THE  MODEL 
7.1  The  Physical  Meaning  of  t 

The  results  given  in  the  previous  sections  show  that  (i)  the  "apparent 
density  of  occupied  trap  sites  at  grain  boundaries  (Nflpp  =  8.5xl021  to 
4.2xl023  nf3)  is  consistent  with  the  range  of  values  that  have  been  pro¬ 
posed  for  the  density  of  trap  sites  at  grain  boundaries  (1019  to  1023  nf3) 
[23];  (ii)  Ngpp  varies  with  pressure,  therefore,  the  parameter  x  is  not 
constant  and  non-equilibrium  partitioning  of  hydrogen  between  lattice  and 
trap  sites  should  be  considered;  and  (ill)  the  exact  value  of  r  cannot 
be  determined  from  Napp  because  Nx  is  not  known. 

Originally,  Bockris  stated  that  t  could  vary  from  0  to  1,  because 
filling  of  the  traps  with  H  was  achieved  by  diffusion  of  hydrogen  from 
the  lattice  into  the  trap  sites.  When  equilibrium  is  not  established, 
then  t  must  be  smaller  than  1  [24].  This  claim  that  x  is  smaller  than 
one  under  non-equilibrium  conditions,  however,  may  not  be  always  true. 

When  hydrogen  enters  the  material  preferentially  through  the  traps,  x 
can  be  expected  to  be  greater  than  one.  Here,  the  establishment  of 

-15- 


'.Vivlv'v 


equilibrium  is  achieved  by  the  transfer  of  hydrogen  from  the  trap  sites 
into  the  lattice.  The  characteristics  of  this  case  may  be  described  as 
follows: 

(i)  Hydrogen  enters  the  material  directly  through  the 

trap  sites  (such  as  prior  austenite  grain  boundaries, 
etc.)  and  then  transfers  into  the  lattice. 

(ii)  Hydrogen  resides  at  the  trap  sites  and  is  prevented 
from  diffusing  into  the  lattice.  As  such,  hydrogen 
concentration  could  build  up  preferentially  at  trap 
sites  and  cause  r  to  become  greater  than  1.  According 
to  Pressouyre  and  Bernstein,  the  binding  enthalpy  to 
the  trap  site,  Hg>60  kJ/mol  dictates  irreversibility  in 
hydrogen  embrittlement  situations  [29], 

(i i 1 )  If  hydrogen  enters  the  material  predominantly  through 
grain  boundaries,  then  the  supersaturability  of  hydro¬ 
gen  should  be  proportional  to  the  hydrogen  supply  rate. 
Thus,  t,  as  a  measure  of  supersaturability,  may  be 
directly  related  to  the  hydrogen  supply  rate. 

From  the  foregoing  discussion,  it  Is  believed  that  the  systems  which 
were  used  to  test  the  model  belong  to  this  case.  This  belief  is  supported 
by  the  following  observations: 

(1)  The  crack  paths  at  low  temperatures  (Region  A)  is 
predominantly  intergranular.  At  higher  temperatures, 
specifically  in  the  transition  region  from  A  to  C, 
the  fracture  mode  is  mixed,  but  with  a  major  part 


by  grain  boundary  separation.  These  grain  boundaries 
are  directly  exposed  to  the  environment  and  are 
expected  to  be  the  preferential  paths  for  hydrogen 
entry  and  diffusion. 

(ii)  The  AISI  4340  steel,  tempered  at  477  K  [6],  contained 

a  high  density  of  Irreversible  traps  in  the  form  of 

high  angle  prior  austenite  grain  boundaries  [30]. 

The  binding  enthalpy  to  the  grain  boundaries  has  been 

estimated  to  be  73  kj/mol,  which  satisfies  the 

requirement  of  binding  enthalpy  for  strong  traps. 

(iii)  The  pressure  dependence  of  t  is  clearly  shown  by 

that  of  N„„.  The  observed  dependence  of  on 
app  r  app 

hydrogen  and  on  hydrogen  sulfide  pressure  (determined 
to  be  0.62  and  0.84  power,  respectively)  compared  well 
with  that  for  surface  reaction  controlled  process 
and  with  that  for  the  transport  rate  of  gases  to  the 
crack  tip,  respectively  [25].  These  correlations 
support  the  earlier  assumption  that  the  non-equilibrium 
parameter  t  is  proportional  to  the  rate  of  hydrogen 
supply  to  the  material. 

Thus,  instead  of  assuming  t  to  be  a  constant  for  a  given  pressure 
is  now  assumed  to  be  proportional  to  the  hydrogen  supply  rate  and 


stVo/T 

Transport 

(12 

esnsP0me*p(-Es/RT) 

Surface  Reaction 

(12 

Woiexp('Ed/2RT) 

Diffusion 

(12 

where  the  g's  are  taken  as  the  proportionality  coefficients.  This  new 
assumption  is  used  to  refine  the  original  model. 

7.2  Refinement  of  the  Original  Model 

Substituting  Eqns.  (12. a)  to  (12. c)  into  Eqns.  (10)  and  (11),  a  new 
set  of  crack  growth  rate  equations  are  obtained  and  are  as  follows: 

1)  For  transport  control 

f  n^t«bfb«(a3/n)Nx(p0/Ti)exp(HB/RT) 


[  l+nt6t6(a3/n)Nx(p0/Ti)exp(HB/RT) 


l+ntet6(a3/n)Nx(p0/T2)exp(HB/RT) 

(P0/Ti> 

2)  For  surface  reaction  control 


fng 6sabfb« ( a3/n )N  Pmexp[ ( HB-Es J/RT] 


l  l+Ssns.S(a3/n)NxP0mexp[(HB-Es)/RT] 


l+n.6.6(a3/n)N  P  mexp[(HR-E  )/RT]J 


P0mexp(-Es/RT) 


l+nd0d6(a3/n)NxPoiexpC(2HB-Ed)/2RT]J 

P0iexp(-Ed/2RT) 


(13. c) 


With  the  Indicated  modifications,  a  new  estimate  of  Hg  of  88  kJ/mol 
Is  obtained  [25].  The  difference  between  this  value  and  the  previous  esti¬ 
mate  of  73  kJ/mol  is  simply  the  activation  energy  for  surface  reaction  (E$) 
or  14.7  kJ/mol  (rounded  to  15  kJ/mol).  Using  this  new  value  of  Hg,  it  is 
seen  that  constant  (though  different)  values  of  Napp  are  obtained  for  hydro¬ 
gen  sulfide  and  for  hydrogen  (with  the  exception  for  the  data  at  13.3  kPa), 
Table  3. 

This  finding  is  gratifying  in  that  the  density  of  trap  site  (Nx)  is 
expected  to  be  a  material  constant  and  the  coefficients  bs,  3t,  ns  and  nt 
are  also  constants.  The  difference  In  Ngpp  for  hydrogen  sulfide  and  hydro¬ 
gen  reflects  differences  in  the  hydrogen  supply  coefficients  (n)  and  the 
coefficients  for  non-equilibrium  partitioning  (&),  and  requires  further 
study.  The  discrepancy  between  Ngpp  for  hydrogen  at  the  lowest  pressure 
(13,3  kPa)  and  those  at  the  other  pressures  may  be  attributed  to  either  an 
underestimation  of  E$  from  the  crack  growth  data  or  the  fact  that  a  transfer 
to  gas-transport  control  had  not  been  considered. 


It  should  be  noted  that  the  modified  equations  produced  no  significant 
improvements  in  predicting  the  crack  growth  response  for  AISI  4340  steel 
in  hydrogen.  Some  improvement  ,  however,  is  obtained  in  the  case  of 
hydrogen  sulfide,  and  is  indicated  by  the  dashed  curves  in  Fig.  6. 

8.  DISCUSSION 

The  model  proposed  in  this  study  has  successfully  combined  the 
achievements  obtained  from  the  micromechanism  st  'dies  [23,  30-32]  and  the 
identification  of  rate  controlling  processes  for  hydrogen  supply  [1,8,9, 
28].  This  unified  "hydrogen  partitioning"  model  successfully  predicts 
the  rapid  reduction  of  crack  growth  rate  associated  with  the  change  in 
fracture  mode  at  high  temperatures.  This  model  may  be  of  help  in  gaining 
Insight  into  the  role  of  microstructural  elements  on  hydrogen  embrittle¬ 
ment. 

As  an  example,  the  effect  of  prior-austenite  grain  size  on  Stage 
II  crack  growth  rate  and  crack  growth  response  may  be  re-examined  in  terms 
of  this  model.  In  the  crack  growth  rate  equations  (13. a)  to  (13. c),  a 
parameter  6  is  Included  and  represents  the  volume  ratio  of  grain  boundar¬ 
ies  to  the  bulk.  If  the  binding  enthalpy  Hg,  the  trap  density  and  the 
other  parameters  are  assumed  to  be  unaffected  by  grain  size.  It  can  be 
seen  readily  that  prior-austenite  grain  size  affects  only  the  temperature 

T  .  at  which  the  crack  growth  rate  attains  its  maximum.  As  shown  in 
max 

Fig.  7,  Ttox  is  observed  to  Increase  with  decreasing  grain  size.  At  room 
temperature,  the  effect  of  grain  size  would  not  be  observed  in  Stage  II 
crack  growth  rate  data  for  grain  sizes  up  to  100  ym,  because  T  would 


be  equal  to  or  higher  than  298K.  When  the  grain  size  becomes  larger  than 


100  urn,  Tmaw  drops  below  298  K  and  the  coarse  grained  material  may  exhibit 
max 

higher  resistance  to  crack  growth.  Fig.  7. 

The  limiting  case  is  that  of  a  single  crystal  where  the  volume  ratio 
6  becomes  zero.  The  crack  growth  rate  is  expected  to  be  much  slower 
because  of  the  absence  of  grain  boundaries.  The  mean  crack  growth  rate 
(da/dt)jj  can  then  be  described  simply  by  Eqn.  (14). 

(da/dt)n  =  a£Q  (14) 

It  should  be  noted  that,  however,  the  contribution  of  quasi -cleavage  was 
neglected  in  deriving  the  crack  growth  equations.  In  the  absence  of 
grain  boundaries,  quasi-cleavage  may  play  an  important  role  in  hydrogen 
embrittlement.  Exclusion  of  its  contribution,  therefore,  may  not  be  rea¬ 
sonable  for  single  crystals  and  for  coarse  grained  materials. 

The  contribution  of  quasi -cleavage  to  hydrogen  embrittlement  can  be 
readily  estimated  from  the  consideration  of  hydrogen-dislocation  inter¬ 
actions.  The  crack  growth  rate  can  then  be  described  as  follows: 

(da/dt)jj  =  (aqfqKq  +  a£f^)Q  (15) 

In  Eqn.  (15),  tcq  and  <£  are  the  hydrogen  distribution  coefficients  for 
quasi -cleavage  and  for  lattice;  respectively,  and  are  given  as  follows: 

5'C  /C£  1 

Kq  =  l+6'Cq/C£  and  Kl  =  l+6'Cq/C£  (16*a) 

Cq/C*  3  T(a3/n)Nxq  exp(Hq/RT)  (16. b) 
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where  NXCJ  is  the  trapping  density  for  quasi-cleavage,  6'  is  the  volume 
fraction  of  quasi-cleavage  sites,  is  the  binding  enthalpy  for  hydrogen 
to  the  quasi -cleavage  planes. 

The  proposed  model  can  be  readily  extended  to  an  examination  of 
the  hydrogen-AISI  4130  steel  system  studied  by  Williams  et  al.  [1-3]*. 
Without  changing  the  binding  enthalpy  Hg,  good  agreement  is  obtained 
between  the  present  model  and  their  experimental  data.  The  comparison 
is  shown  in  Fig.  8,  where  ffa  =  1  and  =  0  are  assumed.  It  should  be 
noted  that  the  form  of  crack  growth  rate  equation  given  by  the  present 
model  (Eqns.  16. a  to  16. c)  and  that  given  by  Williams  and  Nelson  [1-3] 
are  nearly  the  same.  This  similarity  in  form,  however,  does  not  indicate 
agreement,  because  the  physical  bases  of  these  two  models  are  very  differ 
ent.  The  Williams  and  Nelson  model  attributed  the  reduction  of  crack 
growth  rate  at  the  high  temperatures  to  the  reduction  in  hydrogen  supply 
rate.  The  present  model ,  on  the  other  hand,  attributes  this  change  to  a 
partitioning  of  hydrogen  amongst  the  different  microstructural  sites  with 
changes  in  temperature. 

The  effect  of  adsorption  equilibrium,  however,  must  still  be  consid¬ 
ered  as  an  additional  contribution  to  the  present  model.  For  surface 
reaction  controlled  process,  substituting  the  hydrogen  supply  rate,  Q, 
given  by  Williams  and  Nelson  into  the  present  model,  one  can  readily  show 
that  the  crack  growth  equation  at  high  temperatures  reduces  to  the  follow¬ 
ing  form  [25]. 

★ 

Note  that  the  Williams  and  Nelson  data  were  obtained  for  a  fixed  Kt,  and 
do  not  necessarily  correspond  to  Stage  II  crack  growth  over  the  entire 
temperature  range. 


Hr+2aH-2E 

nn(da/dt)n  =  mC  +  RT — (17) 

The  quantity  C  is  approximately  constant  if  the  weak  temperature  dependence 
-k 

of  f^  and  T  ^  is  neglected. 

Let  the  binding  enthalpy  to  the  grain  boundaries,  (Hg),  be  equal  or 
greater  than  58.6  kJ/mol  and  the  apparent  activation  energy  for  surface 
reaction  E$  =  14.7  kJ/mol,  then  the  heat  of  absorption,  aH,  needs  to  be 
equal  to  or  less  than  14.7  kJ/mol  to  conform  to  the  high-temperature  crack 
growth  data.  This  preliminary  result  clearly  shows  that  the  apparent 
heat  of  adsorption  (aH)  is  relatively  small  in  comparison  with  the  binding 
enthalpy  to  the  grain  boundaries.  As  such,  the  contribution  of  adsorption 
equilibrium  to  the  change  in  crack  growth  response  in  Region  C  is  rather 
small. 

The  proposed  model  can  also  account  for  the  change  in  Stage  II  crack 
growth  rate  from  Region  A  to  Region  C  in  the  18N1  maraging  steels.  It 
cannot,  however,  explain  the  sudden  drop  off  in  crack  growth  rates  at  the 
higher  temperatures,  see  Fig.  9.  Mechanisms  involving  two-dimensional 
grain  boundary  surface  phase  transformation,  suggested  by  Gangloff  and 
Wei  [4,5]  and  Chan,  Klier,  and  Wei  [34],  could  provide  a  good  fit  to  the 
experimental  p-T  conditions  for  the  sudden  drop  off  In  Stage  II  crack 
growth  rate.  The  features  of  the  interaction  of  hydrogen  with  particular 
components  of  the  steels  that  account  for  the  specific  responses,  however, 
are  still  unclear.  Further  determination  of  the  relationship  between 
grain  boundary  composition  and  conditions  of  phase  equilibria  involving 
hydrogen  is  essential  to  the  development  of  high  strength  steels  with  high 
resistance  to  hydrogen  assisted  crack  growth  [34], 


9.  SUMMARY 

Results  of  the  current  study  have  provided  an  accurate  characteriza¬ 
tion  of  the  crack  growth  response  for  high  strength  steel  exposed  to 
gaseous  H2  and  H2S.  The  principal  results  from  this  study  are  summarized 
as  follows: 

1.  A  model  which  combines  recent  achievements  from  studies  of  micro- 
mechanisms  of  cracking  and  of  rate  controlling  processes  for  environmen¬ 
tally  assisted  crack  growth  is  proposed. 

2.  The  model  suggests  that  the  main  crack  growth  response  for 
hydrogen-assisted  cracking  is  essentially  determined  by  two  basic  factors: 
(a)  the  rate  of  hydrogen  supplied  to  the  fracture  process  zone,  which  is 
controlled  by  one  of  the  rate  controlling  processes  previously  identified, 
and  (b)  the  partitioning  of  hydrogen  amongst  the  different  microstructural 
elements  or  traps  (principally  between  prior-austenite  grain  boundaries 
and  the  matrix),  which  is  controlled  by  the  hydrogen-trap  interactions 
and  determines  the  micro-crack  growth  rate  for  each  element.  The  differ¬ 
ent  micromechanisms  are  assured  to  operate  in  parallel. 

3.  Quantitative  evaluation  of  the  model  in  relation  to  crack  growth 
response  has  been  carried  out.  By  Including  experimentally  observed 
changes  in  fracture  morphology,  very  good  agreement  between  the  model  pre¬ 
dictions  and  observed  crack  growth  responses  for  AISI  4340  and  4130  steels 


in  H2  and  H2S  is  obtained. 


4.  Partitioning  of  hydrogen  between  prior-austenite  grain  boundaries 
and  matrix  is  shown  to  be  a  more  reasonable  explanation  for  the  observed 
crack  growth  response  at  "high"  temperature,  versus  the  explanations 
offered  by  the  surface  reaction  models. 

5.  Crack  growth  response  is  consistent  with  a  binding  enthalpy  of 
hydrogen  to  the  prior-austenite  grain  boundaries  of  about  90  kJ/mol  (vs. 
about  60  kJ/mol  reported  by  Hirth  for  iron). 

6.  The  contribution  of  prior-austenite  grain  size  to  hydrogen 
embrittlement  was  examined  by  the  model.  The  room  temperature  crack  growth 
rate  is  not  expected  to  be  sensitive  to  grain  size  up  to  100  pm.  For 
grain  size  larger  than  100  pm,  the  room  temperature  crack  growth  rate  may 
decrease  with  increasing  grain  size.  The  fact  that  single  crystals  have 
the  highest  resistance  to  hydrogen  cracking  is  consistent  with  the  model. 


7.  Adsorption  equilibrium  proposed  by  Williams  et  al.  was  considered 
In  relation  to  the  proposed  model  for  surface  reaction  controlled  crack 
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growth.  Preliminary  estimates  Indicated  that  the  apparent  heat  of  adsorp¬ 
tion  may  be  relatively  small  and  Is  less  than  15  kJ/mol .  As  such,  the 
contribution  of  adsorption  equilibrium  is  not  expected  to  be  an  important 
cause  for  the  decrease  in  crack  growth  rates  at  "high"  temperatures. 

8.  Neither  partitioning  of  hydrogen  nor  adsorption  equlllbrlun  can 
account  for  the  observed  "high"  temperature  response  of  18Ni  maraging  steel 
in  hydrogen.  The  two-dimensional  phase  transformation  model  proposed  by 
Chan,  Kller  and  Wei  is  currently  the  only  reasonable  alternative  and  requires 
further  study. 
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endix  I:  Equilibrium  Partition  of  Hydrogen 


Following  Hirth  [23],  the  appropriate  reaction  for  hydrogen  with  a 
single  trap  site  X  is  given  by  Eqn.  (1-1) 


H  +  X  ?  HX 


(1-1) 


The  corresponding  equilibrium  partitioning  of  hydrogen  can  be  expressed 
by  either  the  Boltzmann  [22]  or  the  Fermi -Dirac  Statistics  [23]: 


Hx  »  exp(HB/RT) 


(I-la) 


Hx/(1-Hx)  =  Hl  exp(HB/RT) 


(I-lb) 


where  H  and  H  are  the  concentrations  of  hydrogen  about  the  trap  sites 

A  JC ? 

and  in  the  lattice  (in  atom  fraction)  respectively. 

In  considering  the  partitioning  of  hydrogen  between  grain  boundaries 
and  the  lattice,  the  following  relationships  between  atom  fraction  and 
molal  concentration  is  used: 


CHx  =  <V2VHx 

Ci  *  (No/2NA)Hji  3  E("/a3)/2VHc 


(I-2a) 


( I -2b) 


In  Eqns.  (I-2a)  and  (I-2b),  Nx  is  the  density  of  trap  sites  in  the  grain 
boundaries  (nurrfcers/m3) ;  NQ  is  the  density  of  perfect  lattice  sites  and 
can  be  approximated  by  n/a3,  where  a  is  the  lattice  parameter  of  the  metal 
and  n  is  the  number  of  atoms  per  unit  cell;  and  is  the  Avogadro's  con¬ 
stant.  The  corresponding  equilibrium  partitioning  of  hydrogen  between 
the  grain  boundary  trap  sites  and  the  lattice  (in  molal  concentrations)  can 
be  expressed  as  follows: 


Appendix  II:  Estimation  of  Parameters  from 
Crack  Growth  Data 

The  crack  growth  rate  and  hydrogen  partition  equations  contain  a 
number  of  parameters  that  must  be  prescribed  or  estimated  from  the  experi 
mental  data.  For  the  lattice  parameter  (a),  a  value  of  0.287  nm  for  o- 
iron  is  used.  The  number  of  atoms  per  unit  cell  n  is  equal  to  2.  The 
volume  ratio  between  the  prior-austenite  grain  boundaries  and  the  bulk, 

6,  is  estimated  to  be  10“4,  based  on  an  average  prior-austenite  grain 
size  of  20  ym  for  the  steel  used  In  this  study  and  an  assumed  grain  boun¬ 
dary  thickness  of  1  nm.  The  other  coefficients,  namely  n's,  a‘s  (or  the 
product  net)  and  Hg,  are  estimated  by  matching  the  crack  growth  rate  equa¬ 
tions  to  the  experimental  data  over  specific  temperature  regions.  The 
detailed  estimation  procedures  are  given  in  the  following  subsections. 

Estimation  of  n^ct^,  nsab  and  n^^. 

By  examining  Eqns.  (8a),  (10)  and  (11),  it  can  be  seen  that  the  term 
Cb/C,  ■  (a3/n)Nx  exp(Hg/RT)  is  much  greater  than  one  in  the  low  tempera¬ 
ture  region.  Hence  the  prior-austenite  grain  boundaries  are  essentially 


saturated  with  hydrogen  (<b=l),  and  cracking  is  expected  to  occur  princi¬ 
pally  along  these  boundaries,  with  fb*l.  The  crack  growth  relations, 

Eqn.  (10),  then  reduces  to  the  following  form: 

’  ntab(p/TJs)  Transport  (Il-la) 

(da/dt)jj  a  <  n$otj)pmexp(-Es/RT)  Surface  Reaction  (II-lb) 

n(jabplsexp(-E(j/2RT)  Diffusion  (II-lc) 


These  equations  are  Identical  In  form  to  those  given  by  Wei  [11]. 
Equation  (II-l)  can  be  used  to  determine  the  coefficients  ntab,  nsab 
and  ndab  from  the  experimental  data  on  crack  growth  [25]. 

Estimation  of  nsa£  and  nda£ 

In  the  high  temperature  region,  the  term  Cb/C£  becomes  very  small, 
and  then  approaches  one.  In  other  words,  most  of  the  hydrogen  now 
goes  into  the  matrix  and  cracking  is  expected  to  occur  fully  by  micro- 


void  coalescence  (i.e.,  f ^=1 ) .  In  this  case. 

(da/dt)jj  reduces 

to  the 

following  form: 

VjP/T^) 

Transport 

(II-2  a) 

(da/dt)jj  *  < 

nsa/  exp(-Es/RT) 

Surface  Reaction 

(II -2b) 

ndafcpj£  exP(“E(j/2RT) 

Diffusion 

(H-2c) 

Equatlon  (II-2)  may  be  used  to  estimate  the  coefficients  nsa£  and 

nda£  from  the  aPPr°PrIate  crack  growth  data  [25]. 

Estimation  of  Binding  Enthalpy,  Hg 

It  Is  recognized  that,  at  the  early  stage  of  transition  into  Region  C 
with  increases  in  test  temperature,  the  areal  fraction  of  Intergranular 
(i.e.  grain  boundary)  separation  is  almost  equal  to  1  [16,25].  Also, 
the  rate  of  crack  growth  by  micro-void  coalescence  is  negligibly  small  in 
comparison  with  that  by  grain  boundary  separation  (i.e.,  a£  is  at  least 
3  orders  of  magnitude  below  ab  [25]).  The  contribution  of  micro-void 


coalescence  to  the  crack  growth  rate  (da/dt)jj,  therefore,  may  be  Ignored. 
Equation  (10)  then  can  be  reduced  to  the  following  form: 


(da/dt)n 


< 


’  YbfbT^a3/n'Nx*p/^exp^VRT^ 
VbV6  ( a  Vn  )Nxpmexp[  (H0-Es  )/RT] 

-  ndobfbT6(a3/n)Nxp,sexp[(2HB-Ed)/2RT] 


( 1 1-3  a ) 
( 1 1- 3  b ) 

( 1 1-3  c) 


Equation  ( I 1-3)  can  be  used  to  estimate  the  binding  enthalpy  for  hydrogen 
to  the  prior-austenite  grain  boundaries  (Hg)  from  the  experimental  data 
on  Stage  II  crack  growth  In  the  lower  temperature  portion  of  Region  C. 

For  this  estimate,  Eqn.  (II- 3b)  and  the  data  for  surface  reaction 
controlled  growth  in  an  AISI  4340  steel  at  a  hydrogen  pressure  of  133  kPa 
[7]  are  used.  In  Eqn.  ( II- 3b), all  the  coefficients  except  ffa  and  r 
are  constant  for  a  given  material  and  environment.  Fractographlc  analyses 
have  shown  that  fb  varied  from  1  to  0.8  when  the  temperature  was  Increased 
from  region  A  up  to  the  highest  temperature  tested  in  Region  C  [16,25]. 
Because  of  this  weak  dependence  on  temperature,  fb=l  is  assumed.  The  non- 
equilibrium  coefficient,  t,  may  be  assigned  to  depend  either  strongly  or 
weakly  on  temperature.  In  either  case,  t  can  be  considered  as  a  constant 
for  a  given  pressure  (see  sections  6  and  7).  Finally,  from  a  plot  of 
&n(da/dt)  vs.  1/T  shown  in  Fig.  (II-1)[7],  the  binding  enthalpy  of  hydrogen 
to  the  prior-austenite  grain  boundaries  is  estimated  to  be  approximately 
73  kJ/mol.  This  value  is  higher  than,  but  is  still  comparable  to,  the 
value  of  58.6  kJ/mol  suggested  by  Hirth  [23] . 
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Table  1:  Parameters  Associated  with  Modeling  of  Crack  Growth 
for  AISI  4340  Steel  in  Hydrogen. 


Table  3:  Parameters  Associated  with  the  Modified  Model 
for  Crack  Growth  in  AISI  4340  Steel . 


a.  Hydrogen 


Pressure 

Es 

hB 

Napp  ®snsNx 

kPa 

kJ /mol 

kJ /mol 

( kPa ) _1m“3 

133 

14.7 

88 

3.42  x  1021 

57 

14.2 

88 

3.42  x  1021 

13.3 

8.8 

88 

2.99  x  1020 

b. 

Hydrogen  Sulfide. 

Pressure 

Es 

hB 

Napp  -  etntNx 

kPa 

kJ /mol 

kJ /mol 

K(kPa ) _1m"3 
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